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Recent works have shown the important role Nonlinear Electrodynamics (NLED) can have in two 
crucial questions of Cosmology, concerning particular moments of its evolution for very large and for 
low-curvature regimes, that is for very condensed phase and at the present period of acceleration. 
We present here a toy model of a complete cosmological scenario in which the main factor responsible 
for the geometry is a nonlinear magnetic field which produces a FRW homogeneous and isotropic 
geometry. In this scenario we distinguish four distinct phases: a bouncing period, a radiation era, 
an acceleration era and a re-bouncing. It has already been shown that in NLED a strong magnetic 
field can overcome the inevitability of a singular region typical of linear Maxwell theory; on the 
other extreme situation, that is for very weak magnetic field it can accelerate the expansion. The 
present model goes one step further: after the acceleration phase the universe re-bounces and enter 
in a collapse era. This behavior is a manifestation of the invariance under the dual map of the scale 
factor a{t) — + l/a(t), a consequence of the corresponding inverse symmetry of the electromagnetic 
field {F — > 1/ F, where F = F'^" Fi_iv) of the NLED theory presented here. Such sequence coUapse- 
bouncing-expansion-acceleration-re-bouncing-coUapse constitutes a basic unitary element for the 
structure of the universe that can be repeated indefinitely yielding what we call a Cyclic Magnetic 
Universe. 



I. INTRODUCTION 

In the last years there has been increasing of interest 
on the cosmological effects induced by Nonlinear Elec- 
trodynamics (NLED) [Ij. The main reason for this is re- 
lated to the drastic modification NLED provokes in the 
behavior of the cosmological geometry in respect to two 
of the most important questions of standard cosmology, 
that is, the initial singularity and the acceleration of the 
scale factor. Indeed, NLED provides worthwhile alterna- 
tives to solve these two problems in a unified way, that 
is without invoking different mechanisms for each one of 
them separately. Such economy of hypotheses is certainly 
welcome. The partial analysis of each one of these prob- 
lems was initiated in P, 01 ■ Here we will present a new 
cosmological model, that unifies both descriptions. 

The general form for the dynamics of the electro- 
magnetic field, compatible with covariance and gauge 
conservation principles reduces to L = L[F), where 
F = F^^Ffj^jy. We do not consider here the other invari- 
ant G = F^^F*^ constructed with the dual, since its 
practical importance disappears in cosmological frame- 
work once in our scenario the average of the electric field 
vanishes in a magnetic universe as we shall see in the 
next sections. Thus, the Lagrangian appears as a regular 
function that can be developed as positive or negative 
powers of the invariant F. Positive powers dominate the 
dynamics of the gravitational field in the neighborhood 
of its moment of extremely high curvatures [l| . Negative 
powers control the other extreme, that is, in the case of 
very weak electromagnetic fields 0. In this case as it 
was pointed out previously it modifies the evolution of 



the cosmic geometry for large values of the scale factor, 
inducing the phenomenon of acceleration of the universe. 
The arguments presented in Q make it worth considering 
that only the averaged magnetic field survives in a FRW 
spatially homogeneous and isotropic geometry. Such con- 
figuration of pure averaged magnetic field combined with 
the dynamic equations of General Relativity received the 
generic name of Magnetic Universe 0] . 

The most remarkable property of a Magnetic Universe 
configuration is the fact that from the energy conserva- 
tion law it follows that the dependence on time of the 
magnetic field H{t) is the same irrespective of the spe- 
cific form of the Lagrangian. This property allows us 
to obtain the dependence of the magnetic field on the 
scale factor a(t), without knowing the particular form of 
the Lagrangian L{F). Indeed, as we will show later on, 
from the energy-momentum conservation law it follows 
that H = Hq a~^. This dependence is responsible for the 
property which states that strong magnetic fields dom- 
inates the geometry for small values of the scale factor; 
on the other hand, weak fields determines the evolution 
of the geometry for latter eras when the radius is big 
enough to excite these terms. 

In order to combine both effects, here we will analyze a 
toy model. The symmetric behavior of the magnetic field 
in both extremes ~ that is for very strong and very weak 
regimes - allows the appearance of a repetitive configu- 
ration of the kind exhibited by an eternal cyclic universe. 

Negative power of the field in the Lagrangian of the 
gravitational field was used in Q attempting to explain 
the acceleration of the scale factor of the universe by 
modification of the dynamics of the gravitational field 
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by adding to the Einstein-Hilbert action a term that de- 
pends on negative power of the curvature, that is 

This modification introduced an idea that is worth to be 
generahzed: the dynamics should be invariant with re- 
spect to the inverse symmetry transformation. In other 
words, if X represents the invariant used to construct a 
Lagrangian for a given field, the Action should be invari- 
ant under the map X — > 1/X. Since the Electrodynamics 
is the paradigm of field theory, one should start the exam 
of such a principle into the realm of this theory. In other 
words we will deal here with a new symmetry between 
strong and weak electromagnetic field. In 0, a model 
assuming this idea was presented and its cosmological 
consequences analyzed. In this model, the action for the 
electromagnetic field was modified by the addition of a 
new term, namely 




This action yields an accelerated expansion phase for the 
evolution of the universe, and correctly describes the elec- 
tric field of an isolated charge for a sufficiently small value 
of parameter 7. The acceleration becomes a consequence 
of the properties of this dynamics for the situation in 
which the field is weak. 

In another cosmological context, in the strong regime, 
it has been pointed out in the literature jlj that NLED 
can produces a bouncing, altering another important is- 
sue in Cosmology: the singularity problem. In this arti- 
cle we would like to combine both effects improving the 
action given in Eqn.([T]) to discuss the consequences of 
NLED for both, weak and strong fields. 

It is a well-known fact that under certain assumptions, 
the standard cosmological model unavoidably leads to a 
singular behavior of the curvature invariants in what has 
been termed the Big Bang. This is a highly distressing 
state of affairs, because in the presence of a singularity we 
are obliged to abandon the rational description of Nature. 
It is possible that a complete quantum cosmology could 
describe the state of affairs in a very different and more 
complete way. For the time being, while such complete 
quantum theory is not yet known, one should attempt 
to explore alternatives that are allowed and that provide 
some sort of phenomenological consequences of a more 
profound theory. 

It is tempting then to investigate how NLED can give 
origin to an unified scenario that not only accelerates 
the universe for weak fields (latter cosmological era) but 
that is also capable of avoiding an initial singularity as a 
consequence of its properties in the strong regime. 

Scenarios that avoid an initial singularity have been 
intensely studied over the years. As an example of some 
latest realizations we can mention the pre-big-bang uni- 
verse ITi] and the ekpyrotic universe Q. While these 



models are based on deep modifications on conventional 
physics, that are extremely difficult to be observed, the 
model we present here relies instead on the electromag- 
netic field. The new ingredient that we introduce con- 
cerns the dynamics that is rather different from that of 
Maxwell in distinct regimes. Specifically, the Lagrangian 
we will work with is given by 

LT = a^F^~-F~^ + ^. (2) 

The dimensional constants a, P and jj, are to be deter- 
mined by observation. Thus the complete dynamics of 
electromagnetic and gravitational fields are governed by 
Einstein equations plus Lt- 

We shall see that in Friedmann-Robcrtson- Walker 
(FRW) geometry we can distinguish four typical eras 
which generate a basic unity of the cosmos (BUG) that 
repeat indefinitely. 

The whole cosmological scenario is controlled by the 
energy density p and the pressure p of the magnetic field. 
Each era of the BUG is associated with a specific term 
of the Lagrangian. As we shall see the conservation of 
the energy-momentum tensor implies that the field de- 
pendence on the scale factor yields that the invariant F 
is proportional to a~'^. This dependence is responsible by 
the different dominance of each term of the Lagrangian in 
different phases. The first term a^F^ dominates in very 
early epochs allowing a bouncing to avoid the presence 
of a singularity. Let us call this the bouncing era. The 
second term is the Maxwell linear action which dominates 
in the radiation era. The inverse term fi^ /F dominates 
in the acceleration era. Finally the last term is 
responsible for a re-bouncing. Thus each BUG can be 
described in the following way: 

• The bouncing era: There exists a collapsing phase 
that attains a minimum value for the scale factor 

• The radiation era: after the bouncing, p + 3p 
changes the sign; the universe stops its acceleration 
and start expanding with a < 0; 

• The acceleration era: when the 1/F factor domi- 
nates the universe enters an accelerated regime; 

• The re-bouncing era: when the term 1/_F^ domi- 
nates, the acceleration changes the sign and starts 
a phase in which ii < once more; the scale factor 
attains a maximum and re-bounces 

The universe starts a collapsing phase entering a new 
bouncing era. This unity of four stages, the BUG, con- 
stitutes an eternal cyclic configuration that repeats itself 
indefinitely. 

The plan of the article is as follows. In section II we re- 
view the Tolman process of average in order to conciliate 
the energy distribution of the electromagnetic field with 
a spatially isotropic geometry. Section III presents the 
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notion of the Magnetic Universe and its generic features 
concerning the dynamics of electromagnetic field gener- 
ated by a Lagrangian L = L{F). Section IV presents 
the conditions of bouncing and acceleration of a FRW 
universe in terms of properties to be satisfied by L. In 
section V we introduce the notion of inverse symme- 
try of the electromagnetic field in a cosmological con- 
text. This principle is used to complete the form of 
the Lagrangian that guides the combined dynamics of 
the unique long-range fields yielding a spatially homoge- 
neous and isotropic nonsingular universe. In sections VI 
and VII we present a complete scenario consisting of the 
four eras: a bouncing, an expansion with negative ac- 
celeration, an accelerated phase and a re-bouncing. We 
end with some comments on the form of the scale factor 
and future developments. In appendix we present the 
compatibility of our Lagrangian with standard Coulomb 
law and the modifications induced on causal properties 
of nonlinear electrodynamics. 



II. THE AVERAGE PROCEDURE AND THE 
FLUID REPRESENTATION 

The effects of a nonlinear electromagnetic theory in a 
cosmological setting have been studied in several articles 

Given a generic gauge-independent Lagrangian L — 
L{F), written in terms of the invariant F = Ffj^^F^" 
it follows that the associated energy-momentum tensor, 
defined by 




reduces to 



T., 



7 Sj'^" 



-4:Lp Ffj,"' Fai, — Lg^i,. 



(3) 



(4) 



In the standard cosmological scenario the metric struc- 
ture of space-time is provided by the FLRW geometry. 
For compatibility with the cosmological framework, that 
is, in order that an electromagnetic field can generates 
a homogeneous and isotropic geometry an average pro- 
cedure must be used. We define the volumetric spatial 
average of a quantity X at the time t by 



X 



hm — 

V^Va V 



-5 d x. 



(5) 



where V = j \/—g (Px and Vb is a sufficiently large time- 
dependent three-volume. In this notation, for the elec- 
tromagnetic field to act as a source for the FLRW model 
we need to impose that 



E,Hj = 0, 



(6) 



With these conditions, the energy-momentum tensor of 
the EM field associated to L = L{F) can be written as 
that of a perfect fluid, 



where 



(8) 



(9) 



and Lp = dL/dF. 



III. MAGNETIC UNIVERSE 

A particularly interesting case occurs when only the 
average of the magnetic part does not vanishes and 
E^ = 0. Such situation has been investigated in the cos- 
mological framework yielding what has been called mag- 
netic universe. This should be a real possibility in the 
case of cosmology, since in the early universe the electric 
field is screened by the charged primordial plasma, while 
the magnetic field lines are frozen [3] . In spite of this fact, 
in [2] some attention was devoted to the mathematically 
interesting case in which = a^H^ ^ 0. 

An interesting feature of such magnetic universe comes 
from the fact that it can be associated with a four- 
component non-interacting perfect fluid. Let us give a 
brief proof of the statement that in the cosmological con- 
text the energy-content that follows from this theory can 
be described in terms of a perfect fluid. We work with 
the standard form of the FLRW geometry in Gaussian 
coordinates provided by (we limit the present analysis to 
the Euclidean section) 



ds^ = dt^ - a{tf {dr^ + r^dVi^) 



(10) 



The expansion factor, 9 defined as the divergence of the 
fiuid velocity reduces, in the present case, to the deriva- 
tive of logarithm of the scale factor 



(11) 



The conservation of the energy-momentum tensor pro- 
jected in the direction of the co-moving velocity 
yields 



p+{p + p)G = Q 



(12) 



Using Lagrangian in the case of the magnetic universe 
yields for the density of energy and pressure given in 
equations 



a' F' 



1 



F 



F 



F2 



(13) 



EiEi 



-E gij, HiHj 



(7) 



F' 
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■F 
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where 



F = 2H^ 



(15) 



Substituting these vahies in the conservation law, it fol- 
lows 



a 



= 0. 



(16) 



where Lp = dL / dF. 

The important result that follows from this equation 
is that the dependence on the specific form of the La- 
grangian appears as a multiplicative factor. This prop- 
erty shows that any Lagrangian L{F) yields the same 
dependence of the field on the scale factor irrespective of 
the particular form of the Lagrangian. Indeed, equation 
p6)) yields 



H = Hoa- 



(17) 



This property implies that for each power F'' it is possi- 
ble to associate a specific fluid configuration with density 
of energy pk and pressure pk in such a way that the cor- 
responding equation of state is given by 



Pk 



4fc 



(18) 



We restrict our analysis in the present paper to the theory 
provided by a toy-model described by the Lagrangian 



L1+L2 + L3 + Li 



0" 



a^F^^-F-^^^ 
4 F F^ 



(19) 



where a, (3, p. are parameters characterizing a concrete 
specific model. For latter use we present the correspond- 
ing many-fluid component associated to Lag rangian Lj^ . 
We set for the total density and pressure Pt = Pi and 
Pt = J2 Pi where 



Pi 

P2 



-a" F^ 



1 ^ 1 
— F , Vi — —Pi 

4 ' ^ y 



P3 = ' P3 



Pi 



F2 
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Pi 



Pi- 
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Or, using the dependence of the field on the scale factor 
equation ITTl) . 



Pi 



P2 



P3 



Pi 



- 4a^ H, 

Hol_ 
2 

4 

a 



2H§ 



(21) 



Let us point out a remarkable property of the combined 
system of this NLED generated by Lt and Friedman 
equations of cosmological evolution. A simple look into 
the above expressions for the values of the density of en- 
ergy exhibits what could be a possible difficulty of this 
system in two extreme situations, that is, when F^ and 
1/F'^ terms dominate, since if the radius of the universe 
can attain arbitrary small and/or arbitrary big values, 
then one should face the question regarding the positiv- 
ity of its energy content. However, as we shall show in 
the next sections, the combined system of equations of 
the cosmic metric and the magnetic field described by 
General Relativity and NLED, are such that a beautiful 
conspiracy occurs in such a way that the negative con- 
tributions for the energy density that came from terms 
Li and L4 never overcomes the positive terms that come 
from L2 and L3 . Before arriving at the undesirable values 
where the density of energy could attain negative values, 
the universe bounces ( for very large values of the field) 
and re-bounces (in the other extreme, that is, for very 
small values) to precisely avoid this difficulty. This oc- 
curs at the limit value ps = Prb = 0, as follows from 
equation 



(22) 



We emphasize that this is not an extra condition im- 
posed by hand but a direct consequence of the dynamics 
described by Lt- Indeed, at early stages of the expansion 
phase the dynamics is controlled by the approximation 
Lagrangian Lt ~ ii,2 — Li + L2. Then 



P ■ 



l{l-4a'F). 



Using the conservation law we conclude that the 
density of energy will be always positive since there 
exists a minimum value of the scale factor given by 
^■mim = Sa^H^. A similar conspiracy occurs in the other 
extreme where we approximate Lt ~ ^2,3 = ^2 + ^3, 
which shows that the density remains positive definite, 
since a(t) remains bounded, attaining a maximum in the 
moment the universe makes a re-bounce. These extrema 
occurs precisely at the points where the total density van- 
ishes. Let us now turn to the generic conditions needed 
for the universe to have a bounce and a phase of acceler- 
ated expansion. 



IV. CONDITIONS FOR BOUNCING AND 
ACCELERATION 

1. Acceleration 

From Einstein's equations, the acceleration of the uni- 
verse is related to its matter content by 



3- = -i(p + 3p). 
a 2 



(23) 
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In order to have an accelerated universe, matter must 
satisfy the constraint (p + 3p) < 0. In terms of the quan- 
tities defined in Eqn. ([9]), 

P + 3p = 2{L~4:H'^Lf). (24) 
Hence the constraint (p + 3p) < translates into 

L 



Lf > 



(25) 



It follows that any nonlinear electromagnetic theory that 
satisfies this inequality yields accelerated expansion. In 
our present model it follows that terms L2 and L4 pro- 
duce negative acceleration and Li and L3 yield inflation- 
ary regimes (a > 0). 

For latter uses we write the value of p + 3p for the case 
of Lagrangian Lt ■ 



2. Bouncing 

In order to analyze the conditions for a bouncing it 
is convenient to re-write the equation of acceleration us- 
ing explicitly the expansion factor 0, which is called the 
Raychaudhuri equation: 



1 



1 



(p + 3p) 



(26) 



Thus besides condition ([25)) for the existence of an accel- 
eration a bounce needs further restrictions on a{t). In- 
deed, the existence of a minimum (or a maximum) for the 
scale factor implies that at the bouncing point B the in- 
equality [pB + 'ips) < (or, respectively, (pb + ^pb) > 0) 
must be satisfied. Note that at any extremum (maximum 
or minimum) of the scale factor the density of energy van- 
ishes. This is a direct consequence of the first integral of 
Friedman equation which, in the Euclidean case, reduces 
to equation ((22l) . 



• The inverse symmetry principle: 

The NLED theory should be invariant under the 
inverse map 



F 



For the Lagrangian (|19p . we have chosen the constant 
to be 4/1^. This restricts the number of free parameters 
from three to two, once a direct application of this prin- 
ciple implies that = 16a^/x^. This symmetry induces a 
corresponding one for the geometry. Indeed, the cosmo- 
logical dynamics is invariant under the associated dual 
map 



a{t) 



Ho 1 

a{t) ^ —- 



(27) 



It is precisely this invariance that is at the origin of the 
cyclic property of this cosmological scenario. 

Let us point out that the above map is nothing but 
a conformal transformation. Indeed, in conformal time, 
the geometry takes the form 



ds^ = airif {drf - dr^ - r'^dVl'^ 
Thus making the conformal map 



(28) 



where lo — A/a^, and A = Hq/^/JI. Note that although 
the Lagrangian Lt is not invariant under a conformal 
transformation, the average procedure used to make com- 
patible the dynamical system of the electomagnetic field 
and the Friedman equation is invariant. Indeed, we have 



F = 9 



4^2 

g " Fai_,Fpy - 



V. DUALITY ON THE MAGNETIC UNIVERSE 
AS A CONSEQUENCE OF THE INVERSE 
SYMMETRY 

The cosmological scenario that is presented here deals 
with a cyclic FRW geometry which has a symmetric be- 
havior for small and big values of the scale factor. This 
scenario is possible because the behavior of its energy 
content at high energy is the same as it has in its weak 
regime. This is precisely the case of the magnetic uni- 
verse that we are dealing with here. To obtain a perfect 
symmetric configuration for our model we will impose a 
new dynamical principle: 



VI. A COMPLETE SCENARIO 

Electromagnetic radiation described by a maxwcUian 
distribution has driven the geometry of the universe for 
a period. Let us now analyze the modifications intro- 
duced by the non linear terms in the cosmic scenario. 
The simplest way to do this is to combine the previous 
lagrangian with the dependence of the magnetic field on 
the scale factor. We set 

LT = a^F-~\F^^ + ^^ (29) 

where (3 is related to the other parameters a and p by 
the inverse symmetry principle, as displayed above. 
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A. Potential 

It will be more direct to examine the effects of the 
magnetic universe controlled by the above lagrangian if 
we undertake a qualitative analysis using an analogy with 
classical mechanics. Friedman's equation reduces to the 
set 



+ Y{a) = 



where 



(30) 



(31) 



is a potential that restricts the motion of the localization 
a{t) of the "particle". The constants in V are given by 



B=f, C 
6 



D 



and are positive. 

The dependence of the field as H ~ Hq ja^ implies the 
existence of four distinct epochs, which we will analyze 
now. 

The derivative dL/dF has three zeros, in which p + p 
vanishes. In the case of pure magnetic universe the value 
of F is always positive. We distinguish four eras. 



B. The four eras of the Magnetic Universe 

The dynamics of the universe with matter density 
given by Eqn.(29) can be obtained qualitatively from the 
analysis of Einstein's equations. We distinguish four dis- 
tinct periods according to the dominance of each term 
of the energy density. The early regime (driven by the 
F^ term); the radiation era (where the equation of state 
p = 1/3/9 controls the expansion); the third accelerated 
evolution (where the 1 /F term is the most important one) 
and finally the last era where the 1 /F^ dominates and in 
which the expansion stops, the universe re-bounces and 
enters in a collapse era. 



Using the dependence H 
leads to 



Ho/ a?, on equation pO]) 
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(33) 



We remind the reader that we limit our analysis here to 
the Euclidean section (e = 0). As long as the right-hand 
side of equation ((33|) must not be negative it follows that, 
the scale-factor a{t) cannot be arbitrarily small. Indeed, 
a solution of ([55)1 is given as 



Vlo?). 



(34) 



The radiation period can be achieved from the above 
equation by setting a = 0. As a consequence the av- 
erage strength of the magnetic field H evolves in time 
as 



i/2 = 



1 



2 i2 + 12a2' 



(35) 



Note that at t = the radius of the universe attains a 
minimum value at the bounce: 



(36) 



Therefore, the actual value of depends on Hq, which 
- for given a, fi turns out to be the sole free parameter 
of the model. The energy density p reaches its maximum 
for the value ps — l/QAa^ at the instant t = ts, where 



tB = VUa^. 



(37) 



For smaller values of t the energy density decreases, van- 
ishing at i = 0, while the pressure becomes negative. 
Only for very small times t < ^ Aa^ /k the non-linear 
effects are relevant for cosmological solution of the nor- 
malized scale-factor. Indeed, solution fits the stan- 
dard expression of the Maxwell case at the limit of large 
times. 



2. Radiation era 



1. Bouncing era 

In the strong field limit the value of the scalar of cur- 
vature is small and the volume of the universe attains 
its minimum, the density of energy and the pressure are 
dominated by the terms coming from the quadratic la- 
grangian F2 and is approximated by the forms 



(32) 



The standard Maxwellian term dominates in the inter- 
mediary regime. Due to the dependence on a"^ of the 
field, this phase is defined by » H'^ yielding the 
approximation 



P 
P 



2 

IP 

6 



(38) 



This is the phase dominated by the linear regime of the 
electromagnetic field. Its properties are the same as de- 
scribed in the standard cosmological model. 
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3. The accelerated era: weak field drives the cosmological 
geometry 

When the universe becomes larger, negative powers 
of F dominates and the distribution of energy becomes 
typical of an accelerated universe, that is: 



P 



6 ff2 



(39) 



and we pass from an accelerated regime to a phase in 
which the acceleration is negative. When the field at- 
tains the value Fjib = IGa^fJ.^ the universe changes its 
expansion to a collapse. The scale factor attains its max- 
imum value 

„4 _ ^0 



VII. 



POSITIVITY OF THE DENSITY OF 
ENERGY 



In the intermediate regime between the radiation and 
the acceleration regime the energy content is described 
by the combined form 



or, in terms of the scale factor. 



P ■ 



Hi 



2Hl 



(40) 



For small a it is the ordinary radiation term that dom- 
inates. The \/F term takes over only after a — \/Hq/ ^i, 
and would grows without bound afterwards. In fact, the 
curvature scalar is 



The total density of energy of the BUG is always pos- 
itive definite (see equation [22]) . In the bouncing and in 
the re-bouncing eras it takes the value pB = Prb = 0. At 
these points the density is an extremum. Actually, both 
points are minimum of the density. This is a direct con- 
sequence of equations (fT2l) and (|22p . Indeed, derivative 
of at the bouncing and at the re-bouncing yields 

Ps = ^ p| > 0. 

Thus there must exists another extremum of p which 
should be a maximum. This is indeed the case since 
there exists a value on the domain of the evolution of the 
universe between the two minima such that 



R 



rrifj. 



p-3p 



^0 



showing that one could expect a curvature singularity in 
the future of the universe for which a ^ oo. We shall see, 
however that the presence of the term 1/F'^ changes this 
behavior. 

Using this matter density in Eqn. (P5)l gives 



Pc+Pc = 0. 



At this point we have 



p + Pce^ = 

showing that at this point c the density takes its maxi- 
mum value. 



H§ 1 
3- + -^ ^ 



2 Hi 



0. 



To get a regime of accelerated expansion, we must have 

,2 



3-779 < 0, 



Hi 



which implies that the universe will accelerate for a > Oc, 
with 



Hi 



1/8 



4- Re-Bouncing 

For very big values of the scale factor the density of 
energy can be approximated by 



F 



F2 



(41) 



VIII. THE BEHAVIOR OF THE SCALE 
FACTOR 

Let us pause for a while and describe the form of the 
scale factor as function of time in the four regimes. To 
simplify such description let us separate in three parts: 

Phase A: Bouncing- Radiation 

Phase B: Radiation- Acceleration 

Phase C: Acceleration- Re-bouncing 

characterized respectively by the dynamics controlled by: 
La = L1+L2; Lb = L2+L^\Lc = L3-I-L4. It is straight- 
forward to obtain an analytical expression for each one 
of these periods, which can be analytically continued 
through distinct eras. We obtain for phase A : 



aBR{t) = 



(t-hY 



(42) 
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The inverse symmetric phase C is given by 



a ARB (t) 



2 



[t - + I2a^ 



-1/4 



(43) 



and for the intermediary phase _B, we have: 



aRA{t) 



^2 1 - cos ( Jsn ( (t - ta) , ^ 



1 — cos 



1/4 



where Jsn is the inverse of a first kind elliptic function, or 
the Jacobi function JacobiSN. In order to express the 
analytical continuation of the scalar factor trough the 
different eras it is convenient to re- write equation (|44p in 
the inverse form: 



t - h 



EllipticF 



arccos 



Hl~na^\ ^/2 
Hi + /ia^ 



In the limit jia^ ^ iJg expression (|45l) becomes 

t-t2 



(45) 



2Hl 



Now, in the limit t — ti 12a^ expression aARsit) 
reduces to 



a{t) 



3 " 



-,1/4 



no Big Rip. It describes correctly the radiation era and 
allows for an accelerated phase without introducing any 
extra source. 

The particular form of the dynamics of the magnetic 
field is dictated by the inverse principle, which states that 
the behavior of the field is invariant under the mapping 

F F = ^j^. This reflects on the symmetric behavior 
of the geometry by the dual map a a — Ho/^/JIa. 

The particular form of NLED is based on the princi- 
ple that shows an intimate relation between strong and 
weak field configurations. This inverse-symmetry prin- 
ciple reduces the number of arbitrary parameters of the 
theory and allows for the regular properties of the cosmi- 
cal model. The universe is a cyclic one, having its main 
characteristics synthesized in the following steps: 

• Step 1: The universe contains a collapsing phase 
in which the scalar factor attains a minimum value 

• Step 2: after the bouncing the universe expands 
with a < 0; 

• Step 3: when the 1/ F factor dominates the universe 
enters an accelerated regime; 

• Step 4: when dominates the acceleration 
changes the sign and starts a phase in which a < 
once more, the scale factor attains a maximum and 
re-bounces starting a new collapsing phase; 

• Step 5: the universe repeats the same behavior 
passing steps 1, 2, 3 and 4 again and again, in- 
definitely. 



corresponding to the radiation era, which is continued in 
form URAit) above by identifying ti — t2. 

On the opposite case fxa^ » Hq expression (|45)) goes 

to 



that is 



a{t) 



2 r 



2 fl 

3HI 



{t-hY 



-1/4 



which is continued analytically to the phase C by the 
identification = ti + ttw 



IX. CONCLUSION 

The simplified toy model presented here displays many 
regular properties that should be worth of further investi- 
gation. In particular, it provides a spatially homogeneous 
and isotropic FRW geometry which has no Big Bang and 



Appendix A: Static and spherically symmetric 
electromagnetic solution and the asymptotic regime 

We have made an analysis of the modification of 
electrodynamics in a cosmological context. We are 
not arguing that these effects are more than the re- 
sponse of the universe to local electrodynamics proper- 
ties. Some decades ago Wheeler and Feynman made a 
conjecture that local properties of electrodynamics (e.g. 
the Lienard-Wiechert potential) may just be a conse- 
quence of such cosmic response inducing the elimination 
of advanced fields. However, if one takes these modifica- 
tions as local change of electrodynamics, we should check 
consistency of the theory with conventional electromag- 
netism. We shall restrict ourselves here to the case of the 
static electric field generated by a point charged particle. 
For a general nonlinear Lagrangian L — L{F), the EOM 
for the point charge reduces to 

r^Lp E{r) — const. 

In the case of the Lagrangian given in Equ.® we get 



AE^{r) 



{lQa^E{rf + E{rf - t?E{rf - ff) 



(46) 
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The polynomial in E that follows from this equation 
cannot be solved exactly, but to study the dependence of 
E with r we can plot from (j46|) the function r = r{E). 

It also follows from the plot that E constant for 
r — > oo. 

By taking derivatives of Eqn. (l46|) it can be shown that 
the function E{r) has no extrema ^14;] . Hence, the mod- 
ulus of the electric field decreases monotonically with in- 
creasing r, from an infinite value at the origin to a con- 
stant (nonzero but small) value at infinity. Eqn. (j46|) then 
shows that Eoo = constant ^ 0. This situation is akin 
to that in the theory defined by the action 

It was shown in |^ that the static and spherically 
symmetric solution of this theory does not approach 
Minkowski asymptotically; it tends instead to (anti)-de 
Sitter space-time. We shall see that a similar situation 
occurs in the case of NLED. 

Regarding the behavior of the field for small values of 
r, if we compare the term corresponding to Maxwell's 
case in Eqn. ([351) with the other two, we get that for the 
field to be Maxwell-like there are conditions on the value 
of the free parameters to be fulfilled, thats is: 



« E^ 



In the asymptotic regime, we can set E = Eoo = X and, 
for the energy-momentum tensor one obtains 

TO = Tl = r| = r| = - ^ {X^ + 3^1^ X + 2/?2) 

(49) 

which mimics a A term. 

If we add an extra term in the Lagrangian we could 
eliminate the residual constant field at infinity. In 
the case of Maxwell Electrodynamics such ambiguity of 
choice does not arise due to its linearity. However, for 
non-linear electromagnetic theory a new possibility oc- 
curs which concerns the geometrical structure at infinity. 
This means that for the non-linear electrodynamics the 
fact that at infinity the field is a constant does not im- 
plies that it vanishes. Such property can be translated in 
a formal question, that is, what is the asymptotic regime 
of the geometry of space-time: Minkowski or de Sitter? . 

In classical linear Electrodynamics the answer to that 
question was known and did not pose any ambiguity. No 
longer so if non-linear electromagnetic field is combined 
with the equations of general relativity. The possibility 
of the de Sitter structure must be considered. In theories 
in which a solution distinct from zero for the equation 
Lp — exists, such a question has to be investigated 
combined with cosmology. In a recent paper [l^ a new 
look into this question was considered by the exam of 
a proposed relation of the apparent mass of the graviton 
and the cosmological constant. We will come back to this 
question elsewhere. 



(3^ « E^. 

With the explicit dependence for the field given by 
E{r) = q/r'^, it would be possible to set a value for a in 
agreement with the observation by 

and analogous expressions for fi and /3, where rg is a 
reference value set by the experiment. In the case we fix 
the value of P by the inverse symmetry, the condition on 
f3 reduces to the other two. 



5. Asymptotic regime 

Let US make an extra comment on the above case of 
a point charge particle at spatial infinity. The energy- 
momentum tensor has the form: 

T^. = -Lr^^u - 4.Lf F^a F"i^ (47) 
which in the present case Fqi = Eir) is 

To" Tl ^ {48a^E^ + 2E^ - Q^?E^ - 5/3^) 

T| = T| = - -Ij {\^o?E^ ^ 2E^ + 2n^E^ - /3^8) 



Appendix B: The fundamental state 

A simple look into the equation of motion in NLED 
shows the existence of a very particular solution such 
that its energy distribution is the same as the one in 
the vacuum fundamental state represented by an effective 
cosmological constant. Indeed, the com is given by 

{LpFn-.. = 0. (50) 

Consider the particular solution F = Fq = constant such 
that 

2a^F* - — + fJ-^F ~ 20^ = 0. 

This is the condition that satisfies the equation of mo- 
tion since Lp vanishes at this value Tq. In this state the 
corresponding energy-momentum takes the form 

T^u — A g^in, 

where 

A = -W) = ^Qf3 + ^To-2/3^). 

This property is typical of NLED since there is no possi- 
bility of the linear Maxwell theory having such particular 
solution. 
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Appendix C: Causality 

Most[l^ of our description of the universe is based on 
the behavior of hght in a gravitational field. In a FLRW 
scenario the existence of a horizon inhibiting the com- 
plete interchange of information between arbitrary parts 
of the universe associated with the observation of the 
high degree of isotropy of the CBR generated a causal 
difficulty: different parts of the universe - in the stan- 
dard FLRW geometry could not have enough time to 
homogeneize. The inflationary scenario had its appeal 
precisely for the solution it brought to this problem. 

Non linear theories of Electrodynamics presented a 
completely new look into the causal properties, which 
we will now overview very briefly. 



6. Causal Structure of Nonlinear Electrodynamics 

The main lesson we can extract from the analysis of 
the propagation of the wave front in nonlinear theory is 
contained in a theorem that can be demonstrated using 
Hadamard method to deal with the discontinuity of the 
field and which can be synthesized in a single sentence 
dealing with the modification of the geometry, generat- 
ing an effective metric that controls the properties of the 
space. In order to show this let us make a very short 
summary of it. 



7. The Effective Metric 

We decided to give a very short review of the propa- 
gation of the electromagnetic waves in the NLED, since 
it has very peculiar properties that modify the standard 
description in the linear Maxwellian case. The reader 
interested in more details can consult p^ . 

Historically, the first example of the idea of effective 
metric was presented in 1923 by W. Gordon. In modern 
language, the wave equation for the propagation of light 
in a moving non-dispersive medium, with slowly varying 
refractive index n and 4- velocity u^: 



[d^d^ + {n^~l){u"do.f]F,, 



0. 



Taking the geometrical optics limit, the Hamilton- Jacobi 
equation for light rays can be written as g^'^k^k^ — 
(see 13| for details), where 



(51) 



is the effective metric for this problem. It must be noted 
that only photons in the geometric optics approximation 
move on geodesies of g^^'^' : the particles that compose the 
fluid couple instead to the background Minkowskian met- 
ric. Let us now study in detail the example of nonlinear 
electromagnetism. We start with the action 



where F = F^^Ffj^^, and L is an arbitrary function of 
F. Notice that 7 is the determinant of the background 
metric, which we take in the following to be that of flat 
spacetime, but the same techniques can be applied when 
the background is curved. Varying this action w.r.t. the 
potential W^, related to the field by the expression 



Fn, = W..., - W,,. = W., - W,. 



we obtain the Euler-Lagrange equations of motion 
(EOM) 



{LfF^;. = 0, 



(53) 



where L^? is the derivative Lp = dL/dF. In the particular 
case of a linear dependence of the Lagrangian with the 
invariant F we recover Maxwell's equations of motion. 

As mentioned in the Introduction, we want to study 
the behaviour of perturbations of these EOM around a 
fixed background solution. Instead of using the tradi- 
tional perturbation method, we shall use a more elegant 
method set out by Hadamard. In this method, the prop- 
agation of low-energy photons are studied by following 
the evolution of the wave front, through which the field 
is continuous but its first derivative is not. To be spe- 
cific, let E be the surface of discontinuity defined by the 
equation 

Yi{x^) — constant. 

The discontinuity of a function J through the surface S 
will be represented by [J]e, and its definition is 

— 

The discontinuities of the field and its first derivative are 
given by 



[Fi,u,xh = Uukx, (54) 



where the vector k\ is nothing but the normal to the 
surface S, that is, kx — S^a, and /^i/ represents the dis- 
continuity of the field. 

To set the stage for the nonlinear case, let us first dis- 
cuss the propagation in Maxwell's electrodynamics, for 
which Lff = 0. The EOM then reduces to F.^J" = 0, 
and taking the discontinuity we get 



/^"fc. = 0. 



(55) 



^ L{F) 



(52) 



The other Maxwell equation is given by F*^^'^ = or 
equivalently, 

^^M.;A + F,x■^J. + Fxf,,, = 0. (56) 
The discontinuity of this equation yields 

ftiukx + fuxkfj. + fxfiki, = 0. (57) 
Multiplying this equation by k'^ gives 

U.k'' + Uxk^k^ + fx^k^k, = 0, (58) 
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where k"^ = k^k^^^^ . Using the orthogonahty condition 
from Eqn. ([55)) it foUows that 

f^"'k^ = (59) 

Since the tensor associated to the discontinuity cannot 
vanish (we are assuming that there is a true discontinu- 
ity!) we conclude that the surface of discontinuity is null 
w.r.t. the metric ^'^'^ . That is, 

= 0. (60) 

It follows that k\-fj_k^ — 0, and since the vector of dis- 
continuity is a gradient, 

k^,xk^ = 0. (61) 

This shows that the propagation of discontinuities of the 
electromagnetic field, in the case of Maxwell's equations 
(which are linear), is along the null geodesies of the 
Minkowski background metric. 

Let us apply the same technique to the case of a non- 
linear Lagrangian for the electromagnetic field, given by 
L{F). Taking the discontinuity of the EOM, Eqn. ((55)) . 
we get 

LFf'"'k^ + 2r]LFF F'^'k^^O, (62) 

where we defined the quantity 77 by F"^ fafj = 77. Note 
that contrary to the linear case in which the discontinuity 
tensor /^^ is orthogonal to the propagation vector k^, 
here there is a complicated relation between the vector 
f^'^k^ and quantities dependent on the background field. 
This is the origin of a more involved expression for the 
evolution of the discontinuity vector, as we shall see next. 
Multiplying equation ([SS]) by F^^^ we obtain 

77 fc2 + F^'^Uxk^k^ + F'^-'h^k^k, = 0. (63) 

Now we substitute in this equation the term f^'^k^ from 
Ean. (|62p . and we arrive at the expression 

Lf vk^ - 2 LpFViF^'^k^kx - F^^k^k^), (64) 

which can be written as g^'^k^k^ = 0, where 

g'^'' = LfY' -^LffF^'-F^-'. (65) 

We then conclude that the high-energy photons of a non- 
linear theory of electrodynamics with L — L[F) do not 
propagate on the null cones of the background metric but 
on the null cones of an effective metric, generated by the 
self-interaction of the electromagnetic field. 



8. Causal properties in the fundamental state 

Let us look for the causal structure in the case the elec- 
tromagnetic field rests on its fundamental state. From 
the calculation made in the previous chapter the pho- 
tons propagate in an effective geometry which is given 
by equation (|65p . In the case of the inverse symmetry 
Lagrangian the effective metric tensor takes the form: 



C//) = (-i + f2)5- + ^^-^^.''. (66) 
The fundamental state is the particular solution in which 
F^ = V, 

which corresponds to an energy-momentum tensor equiv- 
alent to a fluid distribution characterized by the condi- 
tion p -\- p = Q and generates a deSitter geometry for the 
background metric g^^ as seen by all forms of matter and 
energy content - as far as we neglect the gravitational in- 
fluence of such remaining matter and energy. However, 
from the above calculation, we conclude that the photons 
do not propagate in such deSitter space but instead in an 
effective metric which is provided by the form: 

This is a very peculiar and interesting situation that can 
be described by the following sentence: 



• The fundamental state of the theory de- 
scribed by the inverse symmetric Lagrangian 
generates a deSitter universe felt by all ex- 
isting matter with one exception: the pho- 
tons which follow geodesies in the above 
anisotropic geometry g'^^f^y 
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